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ABSTRACT 

Recent measurements of the surface magnetic fields of classical T Tauri stars (CTTSs) 
and magnetic cataclysmic variables show that their magnetic fields have a complex 
structure. The magnetic field associated with the octupole moment may dominate 
the magnetic field associated with other moments in some stars, such as the CTTS 
V2129 Oph. Previously, we studied disc accretion onto stars with magnetic fields de- 
scribed by a superposition of aligned or misaligned dipole and quadrupole moments. 
In this paper, we present results of the first simulations of disc accretion onto stars 
with an octupole field. As examples, we consider stars with a superposition of octupole 
and dipole fields of different strengths and investigate matter flow around them, the 
shapes of hot spots on their surfaces, and the light curves produced by their rotation. 
We investigate two possible mechanisms for producing phase shifts in the light curves 
of stars with complex fields: (1) change of the star's intrinsic magnetic field and (2) 
variation of the accretion rate, causing the disc to interact with the magnetic fields 
associated with different stellar magnetic moments. We also discuss the applicability 
of the potential approximation for extrapolation of the surface magnetic field to larger 
distances. 
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1 INTRODUCTION 

Strong magnetic fields have been observed in a number of 
stars, including classical T Tauri stars (CTTSs) (Basri, Marcy & 
Valenti 1992; Johns-Krull 2007), magnetic white dwarfs (e.g., 
Warner 1995; Euchner et al. 2002), and various types of neu- 
tron stars (e.g., Ghosh & Lamb 1978; Ghosh 2007). In most 
stars the structure of the field is unknown. 

It is usually suggested that the magnetic fields of gaseous 
stars are generated and supported by some type of the dy- 
namo mechanism that operates in the stellar interior or in the 
surface layers of the star. Numerical simulations of uniformly 
rotating, fully convective stars show the formation of com- 
plex, non-axisymmetric fields that can be represented only by 
a superposition of different multipoles (e.g., Chabrier & Kiiker 
2006). However, in simulations of differentially rotating stars, 
an additional, ordered component appears (e.g., Dobler, Stix 
& Brandenburg 2006). 
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Measurements of the surface magnetic fields of CTTSs 
using different techniques indicate that they have a complex 
structure (Johns-Krull, Valenti & Koresko 1999; Johns-Krull 
2007). Measurements of the magnetic fields of nearby low- 
mass stars using the Zeeman-Doppler technique show that 
their fields are often complex (Donati & Cameron 1997; Do- 
nati et al. 1999; Jardine et al. 2002). Recent observations 
of two CTTSs have shown that in one star (V2129 Oph) the 
surface magnetic field associated with the octupole moment 
dominates the fields associated with other moments (Donati 
et al. 2007), whereas in the other star (BP Tau) the fields asso- 
ciated with the dipole and octupole moments are both signifi- 
cant and dominate the fields associated with other multipoles 
(Donati et al. 2008). 

Recently, the magnetic field structure of gaseous stars 
with complex fields has been studied analytically. It was found 
that the fraction of the flux through the stellar surface in open 
field lines is smaller in stars with complex fields than in those 
with purely dipolar fields (Gregory et al. 2008; Mohanty & 
Shu 2008). The potential approximation is usually used to ex- 
trapolate the surface magnetic field to larger distances. Gre- 
gory et al. (2006) calculated possible gas flows around stars 
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with magnetic fields constructed from measurements using 
the potential approximation. Recently, we were able to com- 
pute the flow of gas around V2129 Oph and BP Tau in global 
3D MHD simulations in which the magnetic field is not as- 
sumed to be given by a potential but is instead calculated as 
part of the simulation (Long et al. 2009). 

Zeeman tomography of magnetic white dwarfs has 
shown that they also have complex magnetic fields (Euch- 
ner et al. 2002). In some, such as HE 1045-0908, the mag- 
netic field associated with the star's quadrupole moment dom- 
inates whereas the fields associated with its dipole and oc- 
tupole moments are much weaker (Euchner et al. 2005). In 
others, representation of the field requires inclusion of mis- 
aligned dipole, quadrupole, octupole, and other multipoles 
up to n = 4 or n = 5 (Euchner et al. 2006; Beuermann et 
al. 2007). These results are in accord with earlier indications 
of field complexity, such as asymmetric distribution of spots 
on the stellar surface (Meggitt & Wickramasinghe 1989; Pi- 
irola et al. 1987; see Wickramasinghe & Ferrario 2000 for a 
review). 

Neutron stars also have dynamically important mag- 
netic fields. Soon after they have formed, the magnetic fields 
of some may be enhanced by a dynamo mechanism (e.g., 
Thompson & Duncan 1993) and may have a complex struc- 
ture. There is evidence that at least some accretion-powered 
X-ray pulsars have complex magnetic field structures (Eisner 
& Lamb 1976; Gil et al. 2002; Nishimura 2005) or that the 
dipole moment is off-center (e.g., Coburn 2001; see also Rud- 
erman 1991; Chen & Ruderman 1993). Accreting millisecond 
X-ray pulsars have relatively weak dipole magnetic fields (e.g. 
Psaltis & Chakrabarty 1999), but the details of their fields are 
unknown. They have almost sinusoidal light curves that are 
consistent with a dipole magnetic field with a small inclina- 
tion to the spin axis (e.g., Gierliriski & Poutanen 2005; Lamb 
et al. 2008a,b). However, this does not exclude the presence 
of higher-order multipole fields near the star. An unusual phe- 
nomenon — rapid, relatively large shifts in the phase of the 
light curve — has been observed in many millisecond pul- 
sars (e.g., Morgan et al. 2003; Markwardt 2004; Burderi et 
al. 2006; Hartman et al. 2009; Patruno et al. 2009). Several 
mechanisms have been proposed to explain this phenomenon. 

Lamb et al. (2008a,b) have explained several differ- 
ent properties of accreting millisecond pulsars, including the 
phase shifts of their light curves, with the "nearly aligned 
moving spot model". Motion of emitting areas located close to 
the spin axis alters both the shape and the arrival time of the 
pulse in a complex way. These movements could be caused 
by changes in the accretion rate. Earlier 3D MHD simulations 
of accretion onto stars with misaligned dipole magnetic fields 
(Romanova et al. 2003, 2004) have shown that if the dipole 
is nearly aligned with the spin axis, the place where the ac- 
creting matter impacts the stellar surface can move prograde 
or retrograde relative to the surface, which favors the moving 
spot idea. In another model, the magnetic field of the neutron 
star is assumed to change in time (Burderi et al. 2006). Here 
we consider stars with complex magnetic fields and investi- 
gate both of these mechanisms for producing phase shifts. 

In our earlier studies, we were able to perform global 
3D MHD simulations of accretion onto stars with a dipole 
magnetic field (Koldoba et al. 2002; Romanova et al. 2003, 
2004; Kulkarni & Romanova 2005) and onto stars with a com- 
bination of dipole and quadrupole fields (Long, Romanova 




Figure 1. The "cubed sphere" grid of low resolution is shown for 
demonstration. The left panel shows the grid on the surface of the 
inflated cube. Each of the six cube sides has N 2 = 13 2 curvilinear 
Cartesian grid cells. The grid consists of N r = 28 concentric spheres. 
The right panel shows two sectors of the grid. Our simulations use a 
higher grid resolution, up to N r x N 2 = 200 x 51 2 in each sector. 

& Lovelace 2007, 2008). 3D simulations of accretion onto 
stars with multipolar fields are very challenging due to the 
very steep gradients of the magnetic field. In this paper, we 
report on the first simulations of accretion onto stars with 
an octupole field component. Our methods allow us to in- 
vestigate the general case of a magnetic field produced by 
superposition of dipole, quadrupole, and octupole fields ori- 
ented in different directions. However, even the simpler case 
of a magnetic field produced by superposition of misaligned 
dipole and quadrupole fields creates very complex field struc- 
tures (e.g., Long et al. 2008). Consequently, for clarity we 
omit quadrupole fields in the present work, restricting con- 
sideration to the simpler case of a superposition of octupole 
and dipole fields. We investigate matter flows onto stars with 
such fields, the shapes of the resulting hot spots, and the 
light curves produced by these hots spots. We also investigate 
mechanisms for producing phase shifts in the light curves of 
stars with complex magnetic fields. Finally, we discuss the va- 
lidity of the potential approximation for computing the mag- 
netic field outside such stars. 

Section 2 describes the simulation method and the mag- 
netic field configurations. Section 3 presents our results for 
different combinations of dipole and octupole fields. Section 4 
shows examples of the phase shifts that are possible in the 
light curves of stars with complex fields. Section 5 summa- 
rizes and discusses our most important results. 



2 NUMERICAL MODEL AND MAGNETIC FIELD 
CONFIGURATIONS 

2.1 Model 

Our 3D MHD model has been described in a series of papers 
(Koldoba et al. 2002; Romanova et al. 2003; 2004; Kulkarni 
& Romanova 2005; Long et al. 2007, 2008), where disc accre- 
tion onto stars with dipole and more complex magnetic fields 
has been investigated. Here, we briefly summarize different 
aspects of the model. 

2.1.1 Initial Conditions 

We consider a rotating magnetic star surrounded by an ac- 
cretion disc and a corona. The disc is cold and dense, while 
corona is hot and rarefied, and at the reference point ( the 
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inner edge of the disc in the disc plane ), T c = 100T d , 
p c = O.Olpd, where the subscripts 'd' and 'c' denote the disc 
and the corona. The disc and corona are initially in rotational 
hydrodynamic equilibrium, where the sum of the gravita- 
tional, centrifugal, and pressure gradient forces is zero at each 
point in the simulation region. To avoid an initial magnetic 
field discontinuity at the disc-corona boundary, the corona 
is set to rotate with the Keplerian angular velocity at each 
cylindrical radius r. An a-type viscosity is incorporated into 
the code and it operates in the disc only (above some density 
level, p > 0.2), and helps regulate the accretion rate. We use 
a = 0.04 for the octupole-dominated cases, and a = 0.02 
for the dipole-dominated cases. Other parameters of the disc 
such as the density distribution and initial structure are fixed 
in all simulations. 



2.1.2 Boundary Conditions 

At the inner boundary (the surface of the star), most of vari- 
ables A are set to have free boundary conditions, dA/dr — 0. 
The initial magnetic field on the surface of the star is taken 
to be a superposition of misaligned dipole and octupole fields 
(see Bi and B3 in Eqn. 2). As the simulation proceeds, we 
assume that the normal component of the fields remain un- 
changed, i.e., the magnetic field is frozen in the surface of the 
star. We neglect possible changes of the magnetic field struc- 
ture inside the star due to, e.g., dynamo processes, because in 
the most cases the time scale of such variation is much longer 
than the length of the simulation time. 

At the outer boundary, free conditions are taken for all 
variables. In addition, matter is not permitted to flow back 
from the corona into the region. The simulation region is large 
enough [r ma x ~ 36B*), and the disc is massive enough, to 
supply matter for the entire duration of the simulations. 



2.1.3 The "Cubed sphere" Grid 

The 3D MHD equations are solved with a Godunov-type code 
on the "cubed sphere" grid (Koldoba et al. 2002; see also Put- 
man & Lin 2007). The grid consists of N r concentric spheres, 
where each sphere represents an inflated cube. Fig. [l] shows 
that the grid consists of six sectors corresponding to six sides 
of the cube with N x N curvilinear Cartesian grids on each 
side. The whole grid consists of 6 x N r x TV 2 cells. For model- 
ing of octupole fields, a higher radial grid resolution is needed 
near the star compared with the pure dipole cases. We do this 
by choosing the radial size of the grid cells to be 2.5 times 
smaller than the angular size in the region r < (5 — 7)R* 
(see j ]2.1.4| ), while it is equal to the angular size in the outer 
region as in all our previous work. The typical grid used in 
simulations has 6 x N r x N 2 = 6 x 200 x 51 2 grid cells. Sim- 
ulations with higher/lower grid resolutions were performed 
for comparisons. 



2.1.4 Reference Units 

The MHD equations are solved using dimensionless variables 
A. To obtain the_ physical dimensional values A, the dimen- 
sionless values A should be multiplied by the corresponding 
reference units Ao; A — AAo. To choose the reference units, 
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M*(M ) 


0.8 


1 


1.4 


R* 


2R Q 


5000 km 


10 km 


£1* (G) 


10 3 


10 6 


10 9 


Ro (cm) 


4 x 10 11 


1.4 x 10 9 


2.9 x 10 6 


Vo (cm s" 1 ) 


1.6 x 10 7 


3 x 10 8 


8.1 x 10 9 


S2 (s ± ) 


4 x 10~ 5 


0.2 


2.8 X 10 


Po 


1.5 x 10 5 s 
1.8 days 


29 s 


2.2 ms 


B (G) 


43 


4.3 x 10 4 


4.3 x 10 7 


po (g cm" 3 ) 


7 x 1CT 12 


2 x 10~ 8 


2.8 x 10~ 5 


po (dy cm -2 ) 


1.8 x 10 3 


1.8 x 10 9 


1.8 x 10 15 


M °(M yr-i) 


1.8 x 10 19 


1.2 x 10 19 


1.9 x 10 18 


2.8 x 10" 7 


1.9 x 10" 7 


2.9 x 10" 8 


to (g cm 2 s 2 ) 


1.15 x 10 38 


4.9 x 10 36 


4.5 x 10 34 


To (K) 


1.6 x 10 6 


5.6 x 10 s 


3.9 x 10 11 


E (ergs- 1 ) 


4.8 x 10 33 


1.2 x 10 36 


1.2 x 10 38 


T eS , (K) 


4800 


3.2 x 10 5 


2.3 x 10 7 



Table 1. Sample reference units for typical CTTSs, cataclysmic vari- 
ables, and millisecond pulsars. Real dimensional values for variables 
can be obtained by multiplying the dimensionless values of variables 
by these reference units. 



we first choose the stellar mass M* and radius B*. The ref- 
erence units are then chosen as follows: mass Mo = M*, 
distance Ro = B*/0.35, velocity v = (GM /Ro) 1/2 , time 
scale Po = 2-kRo/vq, angular velocity fio = Vo/Ro. The refer- 
ence magnetic field Bo can be obtained by choosing a rea- 
sonable fiducial value for the surface dipole field strength 
Bi*. Then Bo is the fiducial dipole field strength at Bo; 
Bo = Bi*(B*/Bo) 3 . We then define the reference dipole 
moment pi,o = B0B0, quadrupole moment pi2,o = BoRq, 
octupole moment ^3,0 = BqRo, density po = Bq/vq, pres- 
sure po = Po^o> mass accretion rate Mo = po^oB 2 , an- 
gular momentum flux Lo = poVqRq, energy per unit time 
Eq — povoRo (the radiation flux J is also in units of Bo), tem- 
perature To = TZpo/po, where TZ is the gas constant, and the 
effective blackbody temperature T e s,o = (poVo/cr) 1 ^ 4 , where 
a is the Stefan-Boltzmann constant. 

Therefore, the dimensionless variables are r — r/Ro, 

V = V/VQ, t = t/Po, B n = Bn/Bo, Jin = Pn/^n.O (fl = 1,2,3 

for dipole, quadrupole and octupole components) and so on. 
In the subsequent sections, we show dimensionless values for 
all quantities and drop the tildes (~). Our dimensionless sim- 
ulations are applicable to different astrophysical objects with 
different scales. We list the reference values for typical CTTSs, 
cataclysmic variables, and millisecond pulsars in Tab.[l] 



2.2 Magnetic Field Configurations 

When electrical currents outside the star can be neglected, the 
magnetic field there can be described by a magnetic scalar po- 
tential <p(r) and can be represented as a sum of multipoles. In 
this work we consider intrinsic stellar magnetic fields that can 
be described outside the star by a superposition of symmetric 
dipole, quadrupole, and octupole fields. The total magnetic 
field can then be written 

B(r) = B!+B 2 +B 3 , (1) 
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Figure 2. Sketch of the directions of the multipole magnetic mo- 
ments. The rotation axis f2 is in the z direction. The magnetic axes of 
the dipole fix , quadrupole /j.2 and octupole /X3 are inclined at angles 
of ©i, 02 and ©3 from the z— axis. The dipole moment fi± is placed 
in the xz plane. The angles between f2 — /X2, H — fj,3 and xz planes 
are fa and fa respectively. 

where 

Bi = ^llf-^ 

B 2 = g|(5( M ' 2 .f) 2 -l)r-g|(M 2 .r)M 2 

B 3 = 5(||)(/X3-f)[7( / i3-f) 2 -3]f 

-3(||)[5(M3-f) 2 -l]M3. (2) 

Here Bi, B 2 , and B 3 are the magnetic fields produced by 
the symmetric dipole, quadrupole and octupole moments and 
r, fix, (22, and fls are the unit vectors describing the posi- 
tion and the symmetric dipole, quadrupole, and octupole mo- 
ments, respectively. The magnetic moments may be inclined 
at angles 61, 6 2 , and 63 relative to the rotation axis fi. They 
may also be in different meridional planes fi— /x 2 and fl— [13, 
with at azimuthal angles fa and fa relative to the fi — fix 
plane defined by the dipole moment and the rotation axis. 
Fig. [2] illustrates the geometry by an example in which all 
three magnetic moments have different orientations. 

2.3 The Alfven Surface and Magnetospheric Radius 

The magnetospheric radius, r m , is a characteristic radius 
where the inflowing matter is stopped by the magnetosphere. 
Several expressions for this radius have been derived since the 
1970s for stars with dipole fields. The first estimations were 
done for spherically-accreting matter. It was suggested that 
the radially-falling matter is stopped by the magnetic pres- 
sure of the magnetosphere (e.g., Lamb et al. 1973): 

B(r A ) 2 /8tv = p(r A )v{r A ) 2 , (3) 
r A = krf, ^M-^^iGMY 1 ' 7 , (4) 
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Figure 3. The magnetic field distribution on the star with an oc- 
tupole field, fi3 = 0.5, ©3 = 10°, as seen from the equatorial plane 
(left panel), the north pole (middle panel) and the south pole (right 
panel) . 

where r A is the Alfven radius and k is a dimensionless coeffi- 
cient of the order of unity (Eisner & Lamb 1977; Ghosh et al. 
1977). For disc accretion, a similar criterion can be used but 
for stresses. The disk rotates in the azimuthal direction and 
hence main stresses are connected with the azimuthal compo- 
nents of the stress tensor: = [p + pvi] + [B 2 /8tt — B 2 / 47r] 
(here we neglected the viscous stress which is much smaller 
than the matter stress). The motion of the disk is disturbed 
by the rotating magnetosphere, when the matter stresses be- 
come comparable with the magnetic stresses, or p + pv\ — 
B 2 / 8ir — B 2 , / Air . The dominant component of the dipole mag- 
netic field is the poloidal component, and hence B$ << B 
and we obtain the condition for stresses as p + pv 2 , — B 2 /8ir. 

Different criteria were proposed by other authors 
(Cameron & Campbell 1993; Armitage & Clarke 1996; Matt 
& Pudritz 2005). Analysis of these criteria shows that most of 
them give a radius similar to that given by Eqn. [4] but with 
different coefficients k (Bessolaz et al. 2008). A comparison of 
that formula with simulations gives k w 0.5 (Long, Romanova 
& Lovelace 2005), which is very close to the coefficient esti- 
mated by Ghosh & Lamb (1979) . 

For multipolar fields, we use the generalized formula, 
such that the n— th component of the field is B n ~ fi n /r n+2 , 
(see details in Appendix |A2| l, and the magnetospheric radius 

rm = k n r<£l, r { °] n =^ 3 M-^+z{GM)-^+*. (5) 

3 ACCRETION ONTO STARS WITH AN OCTUPOLE AND 
DIPOLE FIELD 

3.1 Pure octupole 

First, we consider accretion onto a star with an octupole mag- 
netic field with p'j, = 0.5. The magnetic axis of the octupole 
field is inclined (misaligned) from the stellar spin axis at a 
small angle, 6 3 = 10°. 

Fig. [3] shows the magnetic field distribution on the stel- 
lar surface. There are two polar regions with strong positive 
(red) and negative (blue) polarities, as seen in the middle and 
right panels. There are also two octupolar belts with opposite 
polarities as seen in the left panel. A meridian line which goes 
from the north to the south magnetic poles will pass through 
the positive northern pole, negative northern belt, positive 
southern belt and negative southern pole with zero magnetic 
field at the magnetic equator. This distribution differs from 
the dipole field which has only positive and negative poles. 

Fig. [4] (top panel) shows 3D views of matter flow to the 
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Figure 4. 3D views of matter flow to a star with an octupolar field 
with fj,^ = 0.5, 63 = 10° at t = 3.5. The disc is shown by a constant 
density surface in green with p = 0. 15 in the top panel. Different den- 
sity levels in the equatorial plane are shown in the bottom panel. The 
colors along the field lines represent different polarities and strengths 
of the field. The thick cyan and orange lines represent the rotation 
and octupole moment axes respectively. 



star at time t = 3.5. The magnetic field near the star has an 
octupolar shape with three sets of closed loops of field lines 
connecting regions of different polarities on the star, which 
we call northern, southern and equatorial sets of loops. Some 
field lines are closed, while others are dragged by the disc, 
wrapped around the rotation axis and inflate into the corona. 
Matter of the disc is lifted above the equatorial set of loops 
and flows into the northern magnetic belt on the star. The 
bottom panel shows the density distribution in the equatorial 
plane, where the disc matter is stopped by the equatorial set 
of loops and forms a low-density gap around the star. 

Fig. [5] shows the density distribution in different planes 
and sample 3D magnetic field lines. Panel (a) shows that at 
t = the magnetic field is a pure octupole field, where the 
three sets of loops of field lines are clearly seen. The equato- 
rial set of loops is expanded to larger distances than the po- 
lar sets, to demonstrate the interaction of the magnetosphere 
with the disc. Panels (b) and (c) show that the disc matter 
is stopped by the magnetosphere, lifted above the disc plane 
and channeled by the equatorial set of loops, and flows into 
the octupolar belts on the stellar surface. Panel (b) shows that 
matter flow is not axisymmetric: more matter flows in from 
the direction towards which the octupole is inclined. Panel (c) 
shows that there is a low-density magnetospheric gap around 
the star. 
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Figure 6. The surface density distribution on the star with an oc- 
tupole magnetic field, ^3 = 0.5, 63 = 10°, as seen from the equato- 
rial plane (left panel), the north pole (middle panel), and the south 
pole (right panel) . 
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Figure 7. The light curves seen by the observer at different inclina- 
tion angles i for the pure octupole case with /13 = 0.5 and ©3 = 10°. 



The bold red line represents the magnetospheric surface, 
where the matter stresses equal the magnetic stresses, 



P 



p + pv 



(6) 



B 2 /8n 

The magnetic field lines stay closed inside the magnetospheric 
surface, where the magnetic stresses dominate. 

The initial octupolar field shown in the panel (a) is a po- 
tential field, that is, a field that has not been disturbed by 
the surrounding plasma. The potential approximation is of- 
ten used in extrapolation of the observed surface fields to 
larger distances (e.g., Jardine et al. 2002; Gregory et al. 2006; 
Donati et al. 2007). Panels (b)-(d) show that the potential 
approximation is sufficiently good inside the Alfven surface, 
where f3 < 1. However, the field strongly departs from the 
potential one at larger distances. 

Fig. [6] shows that the hot spots on the stellar surface rep- 
resent two rings which are located in the planes parallel to the 
magnetic equator. Each ring has a density enhancement on 
one side, where disc accretion is higher due to the inclination 
of the magnetic axis. The middle and right panels show that 
the density enhancements in the northern and southern rings 
are antisymmetric relative to the magnetic equatorial plane. 
That is why the left panel shows only a part of the southern 
ring. 

For the calculation of the light curves, we use the same 
approach as Romanova et al. (2004). The total energy of the 
accreting matter is assumed to be converted into isotropic 
blackbody radiation on the surface of the star. The specific 
intensity of radiation from a position R on the stellar sur- 
face into a solid angle dfi in the direction fc, is I(R, fc) = 
(l/w)F e (R), where F e (R) is the total energy flux of the in- 
flowing matter, 9 = arccos (R ■ fc). Therefore we obtain the 
radiation energy received per unit time in the direction fc, 
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Figure 5. The density distribution in different planes (color background) and 3D magnetic field lines (yellow lines) for an octupole configuration 
with /13 = 0.5, 03 = 10°. Panel (a) shows an xz slice at t = 0; panels (b), (c) and (d) show xz, yz and xy slices at t = 3.5. The red lines show 
surfaces where (3 = 1. The thick cyan and orange lines represent the rotation and octupole moment axes respectively. 

J = r 2 Fobs = J I{R, k) costfdS, where r is the distance be- 
tween the star and the observer, F obs is the observed flux, 
and dS is the surface area element. Simulations show that 
the spots usually "choose" their favorite position on the sur- 
face of the star, which does not vary much with time. That is 
why we fix the spots at some moment of time and rotate the 
star to obtain the light curves. 

Fig. [7] shows the light curves at t = 3.5 viewed at differ- 
ent observer inclination angles i = arccos(S7 ■ fc). They are 
approximately sinusoidal for relatively small inclination an- 
gles, i = 15° and i = 30°, when only the northern ring can 
be seen. For larger i, the southern ring also contributes to the 
light curves and produces a second peak per rotation period. 
The light curve strongly departs from the sinusoidal one at 
larger i (see i = 60°), but two similar peaks per period are 
observed at i = 90°. 



3.2 Strong octupole and weak dipole 

Next, we consider accretion onto a star with a superposition 
of a relatively strong octupole and a much weaker dipole com- 
ponent on the stellar surface: fix = 0.2, /13 = 0.3, 0i = 30°, 
O3 = 20°, cj>3 = 70°. Here, the rotation axis, dipole and oc- 
tupole moments are all misaligned. 

First, we estimate the radius at which the magnetic fields 
of the dipole and octupole components are equal. Using Eqn. 
[2J and suggesting for simplicity that the dipole and octupole 
moments are aligned, we obtain approximate formulae for 
the equatorial and polar magnetic fields: (1) Equatorial field: 
Bi = a*i A 3 . B 3 = (3/2)^3 /r 5 , and (2) Polar field: Si = 
2/ii/r 3 , B3 = 4/i3 /r 5 . By equating Bi and B3 we obtain 
the distances where the dipole and octupole components are 




-200 -150 -100 -50 50 100 150 200 



Figure 8. The magnetic field distribution on the star with a strong oc- 
tupole and a weak dipole component, /ii = 0.2, /13 = 0.3, 0i = 30°, 
03 = 20°, </>3 = 70°, as seen from the equatorial plane (left panel), 
the north pole (middle panel) and the south pole (right panel) . 

equal in the equatorial and polar directions: 

'--(!£)'■ --('£)'• 

Substituting for /ii and we obtain r eq = 1.5 and r po i e = 
1.4. At smaller/larger distances the octupole/dipole field 
dominates. It is evident that near the star, r ~ i?* = 0.35, 
the octupolar component is much stronger than the dipole 
component. That is why in Fig. [8] the magnetic field distribu- 
tion, shows similar features to that in pure octupole case (Fig. 

Fig. [9] shows 3D views of accretion flow onto the star at 
time t = 5. The disc matter comes close to the star and mainly 
interacts with the octupolar component and accretes onto the 
octupolar belts on the stellar surface. External field lines are 
dragged by the disc and inflate. This case shows similar fea- 
tures to that of the pure octupole, because in both cases the 
octupole determines matter flow. 

Fig. [Tot a) shows that initially, at t = 0, the octupole 
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Figure 10. Density distribution in different slices (color background) and 3D magnetic field lines (yellow lines) for the case of an octupole plus 
weak dipole field with parameters fii = 0.2, ^13 = 0.3, 0i = 30°, 63 = 20°, <j>3 = 70°. Panel (a) shows an xz slice at t = 0; panels (b), (c), 
(d) show xz, yz and xy slices at t = 5. The red lines show the Alfven surface, where (3 = 1. The thick cyan, white and orange lines represent the 
rotation, dipole and octupole axes respectively. 



component dominates near the star, while the dipole domi- 
nates at larger distances, such as the inner edge of the disc, 
which is different from the pure octupole case. Panels (b)- (d) 
show that the disc is stopped by the magnetosphere at r m ~ 
0.9 w 2.6-R*. This radius is smaller than r eq « 1.5, and there- 
fore the octupolar field dominates at the disc-magnetosphere 
boundary. The accretion flow is similar to that in the pure oc- 
tupole case (see Fig. [5j, though here the misalignment angle 
= 20° is larger and matter flows more easily above the 
equatorial set of loops. On the other hand, the octupole is in- 
clined more strongly (than in the pure octupole case), and 
part of the funnel stream is located in the equatorial plane, 
due to which the magnetospheric gap is not empty (panel d) . 

Comparison of panel (a) with the other panels shows the 
difference between the initial potential field (panel a) and 
the non-potential fields obtained in the simulations. The mag- 
netic field lines threading the disc and corona inflate, are 
wrapped around the rotation axis and form some kind of a 
magnetic tower (Romanova et al. 2009). The field strongly 
differs from the potential one at r > r m , where the disc and 
corona matter strongly disturb the magnetosphere. The po- 
tential approximation, however, is reasonably good inside the 
magnetospheric radius, r < r m , where (3 < 1 (inside the red 
line in the figure) . 

Fig. [Tl] shows the hot spots on the surface of the star. The 
spots are similar to those in the pure octupole case, but the 
rings have a higher inclination relative to the equatorial plane 
due to the larger misalignment angle. 

Fig.[l2]shows the light curves associated with rotation of 
the star. The shapes of the light curves strongly depart from 
sinusoids, in contrast with the pure octupolar case. However, 
the reason is not the small dipole component, but instead the 
fact that the misalignment angle (03 = 20°) is higher than 



that in pure octupolar case (03 = 10°), and the southern ring 
is seen by the observer even at small inclination angles i. At 
very high i, the observer see both rings in similar proportions, 
and the light curve becomes sinusoidal with two peaks per 
period. 

3.3 Strong dipole and weak octupole 

We now consider the case where the dipole component (/ii = 
2.0) is much stronger than the octupole component (fi 3 = 
0.2). The misalignment angles are 0i = 20°, 3 = 10°, the 
phase angle is ^3 = 180°. 

From Eqn. [7] we find the critical radii where the dipole 
and octupole components are equal: r eq = 0.39 and r po i e = 
0.45. These radii are only slightly larger than the radius of the 
star, R± = 0.35, and hence the dipole determines the dynam- 
ics of the flow in the entire vicinity of the star. 

Fig. [13] shows the magnetic field on the surface of the 
star. The field distribution is complex and different from a 
dipole or octupole field. There are two magnetic poles with 
different polarities near the octupole axis, and two belts of 
different polarities caused by the octupole component. The 
field in the belts, however, is strongly distorted by the dipole 
component. 

Fig. [l4] shows the 3D views of the accretion flow to the 
star at t = 9. The magnetic field looks like a dipole at all 
distances. The matter is stopped by the dipole component of 
the magnetosphere and is channeled to the polar regions in 
two funnel streams which is typical for accretion onto a star 
with a dipole field (Romanova et al. 2003). In the disc plane 
(bottom panel), the matter is stopped by the dipole-like field, 
and no equatorial accretion is observed. 

Fig. [15] (panel a) shows that initially, at t — 0, the mag- 
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Figure 15. Density distribution in different slices (color background) and 3D magnetic field lines (yellow lines) for the case of a strong dipole 
and a very weak octupole field with parameters fi\ = 2.0, (13 = 0.2, ©1 = 20°, 63 = 10°, 03 = 180°. Panel (a) shows an xz slices at t = 0; 
panels (b), (c), (d) show xz, yz and xy slices at t = 9. The red lines show the magnetospheric surface, where f3 = 1. The thick cyan, white and 
orange lines represent the rotation, dipole and octupole axes respectively. 



netic field has a dipole shape in the whole simulation region 
excluding the close vicinity of the star. Panels (b) and (c) show 
that the disc matter is stopped by the dipole-like magneto- 
sphere and accretes towards the poles in two funnel streams. 
Panel (d) shows that a low-density magnetospheric gap forms 
around the star in the disc plane. Comparison of panel (a) 
with panels (b)-(d) shows that the potential approximation 
of the initially dipole magnetic field shown in panel (a) stays 
valid at later times only inside the magnetospheric surface, 
r < r m , where (3 < 1. At larger distances, the field lines are 
dragged by the disc and corona, and the field strongly departs 
from the initially potential dipole field. 

The hot spots are shown in Fig. |16| As we expect, there 
are only two polar spots, which is typical for accretion onto a 
star with a dipole field. No ring-like spots are observed be- 
cause the dipole component dominates almost everywhere 
and matter is not channeled by the octupole-like field. The 
hot spots are located near the northern and southern mag- 
netic poles. 

Fig. ^] shows the light curves associated with the ro- 
tation of the star at t = 9. The shapes of the light curves 
are quite sinusoidal for almost all inclination angles, which 
is also typical for stars with slightly misaligned dipole fields 
(Romanova et al. 2004). At large i, two peaks per period are 
observed, because the spot near the southern magnetic pole 
becomes visible and contributes to the light curves. 

3.4 Dipole and octupole of comparable strength 

We now consider an interesting case where both dipole and 
octupole components disrupt the disc and channel matter. The 
parameters are: ^1 = 1, ^3 = 0.3, 6i = 30°, 9 3 = 20°, 
03 = 70°. For these parameters we find from Eqn. [7j r eq — 



0.67, r po i e — 0.77. However, the inclination of the dipole and 
octupole axes could also influence the result. Fig. ^] shows 
a close view of the matter flow. In the xz plane, the funnel 
streams are first channeled by the dipole-like field. When the 
matter flows close to the star, at about 2i?*, the octupole de- 
termines the flow and converts each funnel stream into three 
small accretion streams between the loops of field lines. In 
the yz plane, the disc is stopped by the dipole component, 
and only a small amount of matter flows along the dipole- 
like field lines and is then governed by the octupole. The xy 
plane shows that matter flow is strongly non-axisymmetric: in 
some directions the disc matter is stopped by the dipole, and 
in other directions matter flows much closer to the star and is 
stopped by the octupolar component. 

Figure [l9] shows the hot spots in this case. Most of the 
matter is governed by the octupolar field and forms ring-like 
spots. Some matter flows into the polar regions located near 
the dipole magnetic poles. 



4 PHASE SHIFTS IN LIGHT CURVES OF STARS WITH 
COMPLEX FIELDS 

In stars with fixed complex magnetic fields, and with a con- 
stant accretion rate, the emitting region may have complex 
shapes, but the location of these hot spots is approximately 
fixed on the stellar surfaces. Rotation of these stars leads to 
a light curve with a pattern that repeats every rotation, and 
no variations like phase shift of the light curves are expected. 
That is what we see in the light curves shown in Fig. [7] |12| 
and |17[ where oscillations occur with the fixed, initial phase. 
However, the pulse from the emitting region may change in a 
complex way and produce the phase shift, if (1) the magnetic 
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Figure 18. Close view of the density distribution in different slices (color background) and magnetic field lines (yellow lines) for the case of a 
strong dipole and octupole magnetic field with parameters fi\ = 1.0, ^3 = 0.3, ©1 = 30°, 63 = 20°, 03 = 70°. From left to right: xz, yz and 
xy slices at t = 3.5. The red line shows the magnetospheric surface, where j3 = 1. 




Figure 9. 3D views of matter flow to a star with a strong octupole and 
weak dipole magnetic field, ^ = 0.2, ^3 = 0.3, 81 = 30°, ©3 = 
20°, <f>3 = 70°, at t = 5. The disc is shown by a constant density 
surface in green with p = 0.2 in the top panel; different density levels 
in the equatorial plane are shown in the bottom panel. The colors 
along the field lines represent different polarities and strengths of the 
field. The thick cyan, white and orange lines represent the rotation, 
dipole and octupole axes respectively. 

field of the star reconstructs, or (2) accretion rate varies with 
time. Below we discuss these two possibilities. 



4.1 Phase shift due to field reconstruction 

The magnetic field of a star may vary with time due to internal 
dynamo processes. Such variations may lead to reconstruction 
of the field, and different multipoles may dominate after the 
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Figure 11. The surface density distribution on the star with a strong 
octupole and weak dipole magnetic field, fii = 0.2, fis = 0.3, ©i = 
30°, ©3 = 20°, 03 = 70°, at t = 5, as seen from the equatorial plane 
(left panel), the north pole (middle panel), and the south pole (right 
panel) . 




Figure 12. The light curves in the case of a strong octupole and 
a weak dipole, ^ = 0.2, p z = 0.3 and ©i = 30°, © 3 = 20°, 
03 = 70°, for different inclination angles i. 
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Figure 13. The surface magnetic field strength of the star with a 
strong dipole and weak octupole magnetic field, pi = 2.0, fj-s = 0.2, 
©i = 20°, ©3 = 10°, 03 = 180°, as seen from the equatorial plane 
(left panel), the north pole (middle panel) and the south pole (right 
panel) . 
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Figure 14. 3D views of matter flow to a star with a strong dipole 
and weak octupole magnetic field, m = 2.0, fi-j, = 0.2, 0i = 20°, 
©3 = 10°, (j>3 = 180°, at t = 9. The disc is shown by a constant 
density surface in green with p = 0.25 in the top panel; different 
density levels in the disc plane are shown in the bottom panel. The 
colors along the field lines represent different polarities and strengths 
of the field. The thick cyan, white and orange lines represent the 
rotation, dipole and octupole axes respectively. 
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Figure 16. The hot spots viewed from different angles for the case 
of m = 2.0, ^3 = 0.2, 6i = 20°, 3 = 10°, 4> 3 = 180°: from the 
equatorial plane (left panel); from the north pole (middle panel); and 
from the south pole (right panel) . The color represents the density. 



reconstruction. Or, the reconstructed field may have the same 
multipoles, but with different orientations of their magnetic 
axes. Here we investigate an example where a predominantly 
octupolar field changes the orientation of its magnetic axis. 
We choose an octupole with fj, 3 =0.3, which in state (a), has 
its magnetic axis at 63 = 20° and <j>-j, = 70°. Later, in state 
(b), it has new orientation angles, O3 = 60° and 3 = 150°. 
In both cases we have a small dipole component: fii — 0.2, 
0i = 30°, which does not change. 

We observe that in both cases the disc interacts with the 
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Figure 17. The light curves in the case of a strong dipole and weak 
octupole, m = 2.0, H3 = O- 2 . ©1 = 20°, 3 = 10°, <fe = 180°, for 
different inclination angles i. 
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Figure 19. The hot spots viewed from different angles for the case 
of Ati = 1.0, ^3 = 0.3, 6i = 30°, 3 = 20°, 03 = 70°: from the 
equatorial plane (left panel); from the north pole (middle panel); and 
from the south pole (right panel) . The color represents the density. 



octupole-like field close to the star (like in §3.1) and the hot 
spots represent two octupolar rings. However, the orientation 
of the rings and their brightness distribution are different due 
to the different axis directions. The most important is that the 
octupole axes have different directions relative to the merid- 
ional plane, which leads to changes in hot spots and that de- 
termines the phase shift. Fig. [20] shows that the light curves in 
states (a) and (b) have different phases with a phase shift of 
A$ = $i, — <!> a — 170°, where the amplitudes are normalized 
to the same values. The shapes of the light curves are differ- 
ent because the angle 3 is different in these states, although 
this difference may not always lead to changes in the pulse 
shapes. 

In most cases it is not well known how fast the dynamo- 
driven magnetic field reconstruction is. Observations of the 
CTTSs type young stars show that the field may possibly vary 
on the scales of months or even weeks (Smirnov et al. 2004; 
Donati et al. 2007, 2008). It is less clear what the time-scale 
of magnetic field variation in neutron stars and white dwarfs 
is. 

We should mention that such phase shifts could occur 
during reconstruction of a field with any combination of mul- 
tipoles, including a purely dipole configuration. 



4.2 Phase shift due to accretion rate variation 

In another example we consider a star with a fixed complex 
field where the different multipoles are oriented at different 
meridional angles, <f> n . Then, at high accretion rates, the disc 
comes closer to the star and interacts with the higher order 
multipoles, while at lower accretion rates, the truncation dis- 
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Figure 20. The phase shift due to the magnetic field reconstruction, 
where predominantly octupolar field in the case of [13 = 0.3, AH = 
0.2, 6i = 30°, changes its orientation. In the state (a) 3 = 20° and 
03 = 70°, in the state (b) 6 3 = 60° and <f> 3 = 150°. 
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Figure 21. The phase shift due to the variation of the accretion rate, 
(a) For low accretion rate (ah = 2), the dipole component stops the 
disc, (b) For high accretion rate (ah = 0.2), the octupole component 
stops the disc. The phase shift occurs during the transition of these 
two states. 



tance of the disc moves away from the star and interacts with 
the lower order multipoles. 

To study this case, we need a slight reinterpretation of 
our reference units defined in |2.1.4| So far, we have used 
the dimensionless parameters /n and /is to control the mag- 
netic moments of the star. However, they can also be used to 
control the accretion rate instead, which is more convenient 
in this subsection; all we need to do is recast the reference 
units into a different form. In order to do that, we note that 
the dimensional dipole moment of the star is ah* — AHAH.o- In 
the previous sections, we kept the reference value /n.o fixed, 
which is why varying ah corresponded to varying the dipole 
moment. Now, however, we keep /h* fixed. It is then con- 
venient to express /j.1,0 in terms of pi* as /h,o = fJ-i*/^!. 
Then, ah,o = Bo-Ro gives B = ah,o/-Ro = aWah.Ro, and 
A*3,o = -Bo-Ro becomes /x3,o = (ii^R^/fXi. The dimensional oc- 
tupole moment then is a«3* = At3A*3,o = (/W^O/Hi-Ro- The 
reference accretion rate is Mo = pqvqRq = BqR^ 2 j\jGM — 

A*?*/ (/ i i-Ro / ' 2 VGM). The dimensional accretion rate Mdim. is 
then given by 

M Mm *(%) ■ (8) 

It is now clear that changing the dimensionless parameter ah 
has the effect of changing the accretion rate Mdim, while 
keeping ah* fixed as mentioned above. To keep a*3* fixed as 
well, we simply have to change a«3 such that the ratio (13*/ /xi* 
is fixed. This allows us to change the accretion rate while 
keeping the stellar magnetic field fixed. 

The physical meaning of this recasting is the following. 
The most natural interpretation of changing ah and a*3 is 
changing the stellar magnetic field, the result of which is a 
decrease in the magnetospheric radius. However, we can also 
decrease the magnetospheric radius by fixing the stellar mag- 
netic field and increasing the accretion rate instead, which is 
exactly what Eqn.[8]shows. The dimensionless parameters ah 



and a*3 are therefore best thought of as controlling the mag- 
netospheric size. 

This brings us the convenience of using ah to represent 
two states of a star with different accretion rates to investigate 
the phase shifts between them. As an example, we take a star 
with a superposition of dipole and octupole fields, where they 
has parameters, 9i = 30°, 9 3 = 20°, (j> 3 = 70° and fixed 
ratio of /is/ah- It is important that the dipole and octupole 
have different phase angles cj> n . 

We choose two states with the above parameters for the 
dipole and octupole but different ah: state (a) with ah = 2; 
and state (b) with ah = 0.2. We observe from simulations 
that in state (a) the disc mainly interacts with the dipole 
component, while in state (b) it mainly interacts with the 
octupole component. Fig. [2TJ shows a phase shift of peaks 
A$ = $b — $ a = 120° between the states (a) and (b), where 
the amplitudes are normalized to the same values. The phase 
shift occurs because the octupole and dipole axes are located 
in different meridional planes and have different angles <f>, 
due to which the matter flows to different places of the star 
at different accretion rates. 

This type of phase shift is expected during periods of 
accretion rate variation. What is the accretion rate "jump" 
between states (a) and (b)? We can give a simple estimate, 
(M) b /(M) a = ((ahW(ah)6) 2 = (2/0.2) 2 = 100. However, 
we should notice that in state (b) octupole strongly domi- 
nates and Eqn. [8] may not be used directly. In addition, the 
transition may occur at smaller values of ah in state (a). So, 
it is clear that this type of phase shift may be produced by 
accretion rate (luminosity) variation. 

Another thing that should be mentioned is that the phase 
shifts would be expected to be accompanied by changes in the 
light curve pulse profile, because the shape and position of the 
hot spots is very different at different field configurations. 
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5 CONCLUSIONS AND DISCUSSIONS 

We performed, for the first time, global 3D MHD simulations 
of accretion onto stars with predominantly octupolar mag- 
netic fields, and for a combination of dipole and octupole 
fields. The calculation become possible due to enhancement 
of the grid resolution near the star. Simulations have shown 
that: 

(i) If the disc interacts predominantly with the octupolar 
field, then matter flows into two octupolar rings on the sur- 
face of the star. In the case of an inclined octupole field, the 
hot spots are inclined and have a brightness amplification on 
one side of the ring. 

(ii) At small misalignment angles of the octupole, 9 = 10°, 
the light curves are quite sinusoidal and "mimic" the shape of 
the light curves of the pure dipole field, in particular at small 
inclination angles i. However, at higher O (e.g., 6 = 20°), 
the light curves strongly depart from sinusoidal even for small 
i. 

(iii) If the dipole component strongly dominates, then mat- 
ter flows in two funnel streams governed by the dipole field, 
and two round spots form in the vicinity of the dipole mag- 
netic poles. If the dipole and octupole have comparable 
strength at the truncation radius, then both components chan- 
nel matter to the star. Usually both octupolar rings and polar 
spots form on the surface of the star. 

(iv) The potential field approximation is valid only in the 
region around the star where the magnetic stresses are higher 
than the matter stresses. At larger distances, the field is 
dragged by the disc and corona and strongly departs from 
being potential. The external magnetic field usually acquires 
a significant azimuthal component and inflates. A number of 
field lines wrap around the rotation axis, forming a magnetic 
tower which may propagate to larger distances due to mag- 
netic force (e.g., Romanova et al. 2009). 

(v) Accretion onto stars with multipolar fields may lead 
to phase shifts in the light curves. This may occur either (a) 
during dynamo-generated internal field reconstruction, or (b) 
during variation of the accretion rate when the complex mag- 
netic field is fixed and has different phase angles <f>„ for each 
multipole component, but the disc interacts with multipoles of 
different orders at different accretion rates. The phase shifts 
would be expected to be accompanied by changes in the light 
curve pulse profile, because the shape and position of the hot 
spots is very different at different field configurations. 

These new challenging 3D simulations helped us under- 
stand accretion onto stars with octupolar fields. Recent mea- 
surements by Donati et al. (2007, 2008) of the two CTTSs 
V2129 Oph and BP Tau have shown that their fields have a 
significant octupolar component. In our next paper (Long et 
al. 2009) we compare our numerical model with observations 
of these stars. 

Magnetic fields in young stars of the CTTS type may have 
a complex structure and may vary with time due to internal 
dynamo processes. The recently measured magnetic fields of 
the accreting T Tauri stars CV Cha and CR Cha (Hussain et 
al. 2009) show a complex structure consisting of a number of 
multipoles. In another CTTS, V2247 Oph, the magnetic field 
is complex and varies on the very short time-scale of a week 
(Donati et al. 2009). Frequent phase shifts are expected in 
this star due to internal field reconstruction, as discussed in 



this paper. In all these stars the phase shifts may also occur 
due to accretion rate variation. 

A phase shift of 0.2 was recently observed in the mil- 
lisecond pulsar SAX J1808.4-3658 on the 14th day of its out- 
burst (Burderi et al. 2006). This phase jump was observed 
only in the fundamental component of the almost sinusoidal 
light curve, and much less so in the 1st harmonic. Burderi et 
al. (2006) suggested that this phenomenon may be connected 
with fast spin-down of the millisecond pulsar or with some 
magnetic field reconstruction during a stage of enhanced ac- 
cretion. We suggest that it raises the possibility that the mag- 
netic field of the millisecond pulsar may be more complex 
than a dipole, and therefore during enhanced accretion the 
disc may interact with deeper layers of the magnetosphere 
where quadrupolar or other components might possibly influ- 
ence the flow. We also suggest that the dipole may be slightly 
off-center (Ruderman 1991; Long et al. 2008), and hence the 
disc matter accretes to both magnetic poles when the accre- 
tion rate is high, or to only one pole if it is low. Superposition 
of dipoles may also be a possible reason. 



ACKNOWLEDGMENTS 

This research was conducted using the NASA High End 
Computing Program computing systems, specifically the 
Columbia, Discover and Pleiades superclusters. The authors 
thank A.V. Koldoba and G.V. Ustyugova for the earlier devel- 
opment of the codes, and A.K. Kulkarni and R.V.E. Lovelace 
for helpful discussions. The research is supported by NSF 
grant AST0709015, and funds of the Fortner Endowed Chair 
at Illinois. Research of MMR is supported by NASA grant 
NNX08AH25G and NSF grant AST-0807129. 



APPENDIX A: THE MAGNETOSPHERIC RADIUS 

The magnetospheric radius r m is the radius where the disc is 
truncated by the magnetosphere. It is also used to depict the 
characteristic size of the magnetosphere. Here we estimate 
the magnetospheric radii for various multipoles. For simplicity 
we assume that they are aligned with the rotation axis. 

Al The Magnetospheric Radius for A Dipole Field 

Here we derive the general formula for the truncation (mag- 
netospheric) radius r m in the cases of the dipole fields. The 
disc is truncated where the magnetic stresses equal the matter 
stresses, B 2 /8tt = p + pv 2 . We assume that the disc is cold, 
p « pv 2 ,, hence B 2 /8ty = pv\. 

We also assume that the system is in a stationary state 
and matter is supplied from the viscous disc with a constant 
accretion rate M — 4ivrhv r p, where r is the distance to the 
star, h = (c s /vK)r is the disc scale height, v r = ac 2 /vK is the 
radial flow velocity, c s is the sound speed, vk = (GM/r) 1 ^ 2 
is the Keplerian velocity, a is the Shakura-Sunyaev viscosity 
parameter. We obtain, 

pvl « pv% = (4 7 ra)- 1 ( l ; fc / Cs ) 3 M(GM) 1/2 r~ 5/2 . (Al) 

For the equatorial component of the dipole field the mag- 
netic field is B — pi/r 3 . Substituting into the equation 
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B 2 /8tv = pv 2 , and using Eqn. Al we obtain, 

r m = farg', fci = (a 2 /2) 1/7 (c s /^) 6/7 , 



(GM) 



-1/7 



(A2) 



The main term, r£? = /4 /7 M~ 2/7 (GM)- 1/7 , exactly coin- 
cides with the main term for spherical accretion (see also 
Bessolaz et al. 2008). However, here the coefficient k\ de- 
pends on a and the ratio c s /vk- We should note that the 
formula for M used in this derivation is applicable only for 
thin discs, and is only approximately valid near the magne- 
tosphere, where the disc matter is accumulated, the disc be- 
comes thicker and has an excess of matter pressure. That is 
why the coefficient ki can be considered to be approximate 
and can be substituted with some number. Comparison of nu- 
merical results for r m with Eqn. |A2| made in Long et al. (2005) 
shows that for their set of simulations, fe a « 0.5 (see more 
comparisons in Bessolaz et al. 2008). 

A2 The Magnetospheric Radius for A Multipolar Field 

Next, we derive the magnetospheric radius for multipolar 
field. For simplicity, we assume that the main axis of the mul- 
tipolar field is aligned with the rotation axis, and therefore 
the n— th multipolar component of the magnetic field can be 
presented in the form B n ~ fj, n /r n+2 . Similarly, we obtain, 
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This formula can be applied in the cases of aligned multipoles, 
or can be used for estimates of the magnetospheric radius in 
general, misaligned cases. 

We should note that the B z equals in the disc plane for 
aligned fields of 2n-th order multipoles, such as quadrupole 
(ref. Long et al. 2007), and matter could flow directly to the 
star in the disc plane. So r m ,„ does not reflect where the in- 
flowing matter stops, but only the size of magnetosphere. 
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